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ABSTRACT 

The pressure and phase distribution of sound in a fast fluid 
medium underlying a tapered fluid medium was modeled using a Green's 
function approach. The model predicted well defined beams in the 
bottom, with apparent interference effects strongly evident for cases 
at a distance from the wedge interface. The effects of attenuation on 
the patterns were studied. A simple expression for determining the 
beam depression angle as a function the wedge angle, the sound speed 
ratio, and the density ratio was derived. Comparisons with laboratory 


measurements were made. 
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I. INTRODUCTION 


The problem of radiation behavior in a wedge-shaped medium 
overlying a bottom possesing a greater speed of sound has been 
investigated both optically and acoustically. In 1971, Tien and 
Martin (Ref. 1) deposited a thin, tapered, dielectric film ona 
substrate having a higher refractive index, and observed the behavior 
of a laser beam coupled into the film. Perfect reflection at the 
film-substrate interface was observed until the changing angle of 
incidence decreased below the critical angle for continued mode 
propagation. For angles of incidence below critical, the light was 
converted into radiation modes in the substrate. 

A concurrent analysis of the acoustic problem was done by Kuznetzov 
(Ref. 2). Kuznetzov developed a theory, based on normal~modes, that 
predicted the following for sound propagating within a wedge-shaped 
medium overlying a half~space with a higher sound speed: 

l. As the sound approached the vertex of the wedge, it would 
continue to be totally reflected until the angle of incidence decreased 
temene limiting angle of total reflection. 

2. Sound incident at less than the limiting angle of total 
reflection would be totally refracted into the underlying half-space. 
This refraction occurs within the region from the wedge apex to the 


point where the limiting angle was attained. 





3. Within the half-space, the acoustic energy would be propagated 
as a well defined beam with the maximum pressure amplitude Occurring at 
an angle of depression from the wedge interface within the range B to 28, 
where 8 is the wedge angle. 

Kuznetzov also performed a series of experiments that supported his 
EMCOrY. 

Subsequent optical experiments by Tien, Smolinsky, and Martin (Ref. 3) 
demonstrated the formation of a well-defined beam within the substrate. 
The accompanying theoretical discussion predicted conflicting behavior. 
Ray-optics predicted refraction into the substrate after the cutoff dis- 
tance (the point where refraction begins) and wave theory predicted com- 
plete refraction into the substrate before the cutoff point. Sigelman, 
et. al. (Ref. 4) both predicted and observed a well columnated beam in a 
water-aluminum wedge system, but the observed beam was much broader than 
predicted. Also, the pressure maxima for the columnated beam occurred at 
approximately 11° for a wedge angle of 1.3°, well outside the limit pre- 
dicted be Kuznetzov. 

The primary objective of this theoretical development is to produce 
a simple, fast computer model to predict the pressure and phase distri- 
bution of the columnated sound beam in the fast substrate. The pressure 
and phase distribution along the wedge interface calculated by Kawamura 
(Ref. 5) will be used as input to the substrate model, and will specifi- 
cally be used as the boundary conditions for a Green's function analysis. 
To remove explicit frequency dependence, all relevant distances will be 
normalized to the distance measured from the apex along the wedge inter- 


face at which the critical reflection angle is first exceeded; this dump 





distance is referred to as the dump distance X. Integration over the 
wedge interface will be truncated at distances sufficient to encompass 
only selected modes of interest. Results will be analyzed for trends 


and compared with existing measured data. 
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Pressure and phase distributions along the bottom of a wedge- 
shaped medium are obtained using an adaptation of the method of images 
as prescribed in Ref. 5. The source may be set at either finite or 
infinite distance, but must be maintained sufficiently far from the 
wedge boundaries to ensure that impinging waves are essentially locally 
planar. The following theoretical discussion uses the above pressure 
and phase distribution as initial conditions, and calculates the 
pressure and phase distribution in the fast bottom by utilizing the 
Green's function (Ref. 6). 

Figure II-1 illustrates the geometry of the problem. The horizontal 
X-axis and the vertical z-axis have their point of origin at the wedge 
apex, with the positive x-axis forming the interface between the wedge 
and fast bottom. The wedge angle is designated 8 and the angle 
Beccieying the sound source location is ¢. X is the first mode 


dumping distance, and is a distance determined as follows: 


o. SEaAreGOs (Co 1/7 Cin) (II-la) 
DL = a / { 4-f-sin(0.] } (II-1) 
Koo DA / sin(8) (II-1c) 
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where 


C; = sound speed in the wedge 
Co = sound speed in the bottom 
f = frequency 


and D. = lowest mode cut-off depth. 
To make the problem frequency independent, all distances were 


normalized to X by defining 


ee (II-2a) 
emes (x Vex (II-2b) 
and fo 32 (II-2c) 


The variable x” is identical to x, and &* is identical to —. The 
prime designation is required for dummy variables of integration in 
the Green's function development to follow. 


Pressure along the wedge bottom is then given by 


p(—,n=0,t) = P(E) exp fifut + L(E)}} (II-3) 


where P(&) is the pressure amplitude along the wedge bottom and L(&) 
is the pnase along the wedge bottom. Pressure in the bottom, away 
from the boundary, is then given by Ref, 7, 


X az’ (II-4) 
7 =() 


p(r,t) 
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exp[iut] = ( P(E’) a 


exp [ iwe | R(r) 
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where G is the appropriate Green's function. The z° coordinate is 
. : 7; > > -~ 
identical to the z coordinate and r = (x? + y2y 7, The appropriate 


Green's function, for a line source in cylindrical geometry is 


—_ 


G(r,r°) = inn?) (ko |r - 271) (II-5) 


4 : Bie ReaD ante 
where Ka 1s the propagation vector in the bottom and r = [ (x erect az) ie 


Differentiation of the Green's function yields (Ref. 7) 


ko (Zz = 2%) 
Genera) emer Pay 2 
ee a (ky |r ne 1) eas Clie) 


=> re 2) & 
where ko|x = | >> 27 and He ey H. are zero and first order Bessel 
Functions of the third kind (Hankel functions) respectively. Using 
: , —_ =e 
equations (II-2a,b,c) at z~ = 0 the expression lx = may be 


rewritten as 
ie => ano 2)% 
\r- 2 | = X [ (e - €°)* +n (TT=7) 


: G : 
Substituting (II-7) into equation (II-6), and evaluating poerae Ze) 












gz 
yields 
t it Mes 
eu T exp(=*} exp{-ikox[ (E-£7) 2 + n?| “}konX 
— — II-8 
< 4 
OZ a ela ae : nz" 


Combining equations (II-8) and (11-4) yields a time independent 


Se 
expression for amplitude and phase, R(r) 


1 
ww ' 2k5X 
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P(E“ )expli[L(E~) -kox [(E-67) 7407] 77D ab" (II-9) 
a 273/4 
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Equation (II-9) is the basic equation used for the calculation of 





Phase and amplitude in the bottom. 
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IftI. DEVELOPMENT OF WDGBTM 


Program WDGBTM was designed to perform four main functions. 

First, the pressure and phase along the interface of a wedge and 
underlying fast media may be calculated in two ways, for both 
infinitely and finitely distant sound sources in accordance with 

Ref. 5. Modifications were made to the code of Ref. to facilitate 
subsequent in-bottom calculations, and to enhance user flexibility. 
These changes did not affectresults. Then using the first code as 
input, additional code enables both the computation of pressure and 
phase along a line perpendicular to the wedge bottom, and as a function 
of radius and angle measured from the first mode dumping point, X. 

Modes of program execution may be chosen by literal input directives 
combined with appropriate numerical input. Appendix A. illustrates an 
input set that utilizes all four options. Output consists of both 
printed listings and graphics. 

Program versions exist for execution on both a CDC 6500 and an IBM 
360. The CDC 6500, with a 60 bit word, does not require double 
precision for accuracy equivalent to the IBM 360 output. The CDC 
version of WDGBTM has simpler varian control routines, and separate 
output files for each pressure/phase variation. This enables optional 
Suppression of listings. Execution time on the CDC 6500 is roughly ten 
times faster than on the IBM 360. For average input sets, times are 150 


and 1500 cp seconds respectively. 
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A. EQUATION PREPARATION 
leben Domain 
Integrations along lines perpendicular to the wedge bottom 
were designated depth domain calculations. The program was designed to 
allow specification of the field point, i.e., specification of where 
along the media interface line the perpendicular probe line could be 
dropped. The field point is not confined to lines directly beneath the 
wedge, but may extend to points beyond the wedge apex. These points 
would be values of negative * in Figure II-l. The beginning and ending 
points of the section of the perpendicular to be examined, and the 
number of field points used in the interval were also to be specified. 
Equation (II-9) was split into real and imaginary components before 
coding. Integration order selected was over the entire E~domain for 
each.) point. 
2. Angle Domain 
Examination of FigureII-l1 verifies the following geometric 


relationships. 


Mm 
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O cos0,, (III-la) 


and 
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where € is the distance from X to the &~ coordinate of the field point 


of interest and © is the radial distance from X to the field point. 
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Substitution of equations (III-la) and (III-1lb) into equation (II-9) 





yields eee 
Pa (G76) i int \ (2k X 1% 
= i exp ( - (III-2) 


P(E“ expt if L(E*)-kox[ (E-67)?+07-e7] 1/7] jae’ 
Equation (III-2) is the basis for angle domain calculations. The 
specification of o originating at X was chosen to conform to existing 
experimental data (Reference 8). The length of the radius, the 
beginning and ending of the angular sweep, and the angular increment 
comprise the input specifications. As in the depth domain option, 
equation (III-2) was split into real and imaginary components before 
coding. 
Se Design Restrictions 

Theoretical development of program WDGBTM assumes the sound 
source is sufficiently distant to ensure essentially locally planar 
waves at the interface and the in-bottom measurements satisfy the 
condition Ke |r-z >>2T. 

4. Program Results 

An initial set of WDGBTM runs were made in the depth domain for 
a perpendicular directly down from the wedge apex. The source was 
infinitely distant and set at an angle of B/2. This limited modes within 
the wedge to odd numbered (symmetric) modes only. Integration over the 
wedge bottom was truncated at § equal four, so that only modes one and 
three were considered as contributing significantly to the beam pattern 
in the bottom. Beam patterns in the bottom were calculated over the range 


inmecueal zero to 1.3X. 
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The Initial set was comprised of seventy-five runs. In addition 
to the above criteria, the runs were generated by all possible combi- 
nations of the following specifications. The ration 0/0, was 0.95, 
0.90, 0.80, 0.70 or 0.50. The sound speed ratio, C,/C. was set at 0.95, 
0.90, or 0.85. The wedge angle 8 was varied through 10, 6.92, 5, 2 and 
Il degree(s). 

From the results of the initial data set, the depression angle and 
beam width was calculated for the beam pattern in the bottom. Results 
are summarized in Appendix B. These results were then used to find a 
simple analytical formula to predict the depression angle as discussed 
maeseceion IV. 

Other runs were made for perpendiculars at x =-9X from the wedge 
apex. These runs displayed apparent strong interference phenomona and 
are discussed in detail later. 

Angle domain executions were made to compare with experimental data. 


Results are discussed in Section VII. 


BemomnionrrViTyY TO INTEGRATION STEP SIZE 

It was desirable to keep both the integration step size and the field 
point increment size as large as possible in order to minimize total 
computer resources required. Increment of 0.01 times the dump distance 
(X) proved to be optimum for both integration increments along the 
wedge interface and field point steps in the bottom. 

To illustrate the adequacy of 0.01X increment, a series of runs were 
made in accordance with Table III-1l. An overlay of wedge interface 


pressure amplitude for AX = 0.01X and AX = 0.001X demonstrated no dis- 


tinguishable difference. Likewise, depth domain plots (a, b, c, d of 


Lo 





Table CII - l) 


Sensitivity to Integration Step Size (Ax) and Field point interval 


size (An) 
me = 0.002 Ax = 0.0001 
An=0.01 An=0.001 An=0.01 An=0.001 

(a) (b) (Cc) (d) 

Run Nhamp max on 
a 0.240 £3.50 
b 0.244 aah 
on 0.240 ets Sie 
da 0.243 i506 


oe7, 





Table IIfi-1) could be 1dentically superimposed. Use of an integration 
step size of 0.001X resulted in a calculated depression angle within 

one percent of the depression angle for the 0.01X run. Utilization of 

a step size of greater than 0.01X is not recommended. Within this 
limitation the integrations can be carried out over 10X along the inter- 
face for a maximum field interval of 10X in the bottom. The 10X re- 
striction is imposed by program array size. (For the 0*s specified, 0.5 
degrees in the angle domain was slightly finer than the 0.01X field 
point step size in the depth domain and was the field point increment of 


choice). 
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IV. BEAM DEPRESSION ANGLE ANALYSIS 


A simple relationship was sought between the physical parameters 
associated with the wedge system and the resultant beam. 


Assume 

eeeca iva) 
where a is the beam depression angle and 8 is the wedge angle. For 
each 94/0, and c,/C, 1 GoaVvs 8 was graphed on log-log paper, and the 
Slope (n) was measured. Results are listed in Table 1 of Appendix C, 
meeequals 0.329. 


Then Onc Was found as a function of c,/e Dhe critical angle, § 


25 ee 


Menon explicit function of c,/C. 


90 = arccos (C,/C,) 9c in radians Civ 2) 


It was then assumed 


9 a ee 9 M 


nC Cc (IvV-3) 


For any constant 8, M was derived as the slope of a log-log graph of 


vs 9 Graphs were analysed for 8 equal to 10, 6.92, 5, 2 and l 


Joe Ce 


degree, and p,/79. ratios of 0.95, 0.90, 0.85, 0.70 and 0.50. Results 


are summarized in Table (2) of Appendix C, “% = 0.772. 


ey) 








The same process was repeated for ascertaining dependence on 


P,/P. - Log-log graphs were prepared by holding 8 and Fa constant 


using the assumed relation 
0.329 Os7 72 


eee Cr 


q 


Results are summarized in Table 3 of Appendix Cg = -0.254. 

The proposed formula is now dependent on all physical properties, 
assuming the source angle is maintained at one-half 8. To remove the 
proportionality sign, write 

OG 2329 O.7 72 ~0.254 


eas B Bn (0/95) (IV-5) 


and solve for K. From Table 4, Appendix C, K = 17.22/radian. 
Therefore 
Oia! OL77 2 -0.254 
= lV- 
B nc Dyes OO (0 ,/P >) ( 6) 
where 8 is in degrees ~d Soin radians. 
One calculated from equation (IV-6) and the difference between J ne 


respectively) are also 


and 85 from program WDGBTM (labeled AG, and Oe 


tabulated on Table 4, Appendix C. Ag, varies from a minimum of 0.00 to 
a maximum of 1.25, with the average of the absolute values of AG , equal 
to 0.36. Errors ranged from zero percent to 5.0 percent. Thus oe gives 


a reasonable approximation of the beam depression angle. 


20 





aesimular f£Lormula 


i 1/3 OST Ts =0.27 agp) 
Opa ee nn oa (0 ,/0 5) 
was developed independently by A. Coppens (Reference 6) The minimum 
deviation of 0 5, trom J, 1S 0.01 and the maximum deviation is 1.56. 
The average of the absolute value of the deviation is 0.36. Both 


formulas are best for 01/9, close to 1.0, and less accurate for 


01/05 equal to 0.5. 


Zi 
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V. THE EFFECT OF ATTENUATION 


As previewed in Section III-A.4 the pressure distribution as a 
function of N at x = -9.0X demonstrated apparent interference effects 
(see figure V-1). As integration over the wedge interface to 4X 
encompasses dumping by the first and third mode, and the third mode 
dumps about the point 3X vice X for the first mode, it is logical to 
expect attenuation effects to reduce the contribution of the more 
distant mode. The initial Green's function approach, equation (II-9) 
did not consider attenuation of the acoustic wave in the bottom 


medium. An attenuation term 


exp {-a ioe +t naleaee (V-1) 


was incorporated into equation (II-9) to yield 


ait = 1/4 exe(2) (2824) (y=2) 





P(E’) exp{—af (E-E*) “en? ) 2/7 texpla[ (67) -kax [(B-6 7) 74? Jo 2) Fab" 
5 eae) a 


An appropriate code modification to the depth domain routine was made. 
The attenuation coefficient % was an input variable to the program. 

Two sets of comparison runs were made. The first set consisted of 
an unattenuated control (®= 0) and an attenuated run with © equal 1.0, 
with the pressure amplitude measured down from the wedge apex. The 
second set consisted of an unattenuated control plus runs with @ set to 


1.0, 0.5, 0.1, and 0.01 and the pressure measured at * = -9.0 from the 
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apex. All other input variables for both sets were held constant and 
are specified in Table (V-1). Integration over the wedge bottom was 
from the apex to four times the dump distance; this was to insure that 
contributions from the third mode (dumping about the point x= 3X) were 
included. 

The results of set one are illustrated in Figure (V-2). Note that 
the second maxima is severely depressed in the attenuated case. The 
depression angle has been shifted slightly from 11.86 degrees for the 
unattenuated data to 11.31 degrees with © set to 1.0. 

Figure V-3 is a plot of the results from set two. All curves have 
been normalized to the initial slope for = 1.0. The second and third 
maxima decrease with increasing %, and relatively more energy is 
contained in the first beam. The depression angle of the first beam 
decreases slightly with increasing ®@ . These trends are summarized in 
Table (V-2). 

Figures (V-4a) and (V-4b) illustrate the effects of attenuation on 
different contributing modes. For Figure (V-4a), the source was placed 
at 1/3 8. With this source geometry, contributions from mode three are 
suppressed, and only the first and second mode contributed to the beam 
pattern in the bottom. For Figure (V-4b), the source was placed at 58 
and modes one and three propagate into the bottom medium. Both cases 
were run for @ = 1.0, 0.5, 0.1 and 0.0. As could be expected, 
attenuation causes greater suppression of third mode interference 


contributions than of second mode contributions. 
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Table (V-1) 


Input Specifications for Attenuation Runs 


= IE se = P 
Cc, 500.0 Ps lO 
C. = 1578.947 Pe. =k O52652 
Cie = 0.95 PA, = 0.95 
= 5.0° he—=) 5000 0rrz 
xX, = 4.0 Source distance = infinite 
§ =2.5 
(Oe O-r- OO Oe LO 5 bo 


Refer to Figure II-l for variable meaning. 
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VI. ANALYSES OF FAR-FIELD INTERFERENCE PATTERNS 


The results of the attenuation investigation were supportive of an 
interference explanation for the multiple beams. Subsequently, two 
analytical approaches were tried to predict the location of interference 
maxima and minima under unattenuated far-field conditions. 

Both the Differential Change and Multiple Lambda analyses assume 
very simple processes. Both assume point sources of sound at the wedge 
interface located at X and 3X. Reallistically the sound sources behave 
as phased arrays propagating assymetrically about the dumping points. 
Phase differences between the two dumping points are non-zero. The 
Differential Change formula becomes invalid as approaches Zero. 

A. DIFFERENTIAL CHANGE APPROACH 
Based on the simple geometric considerations illustrated in Fig (VI-l), 
the following nondimensional formula for the change in interference 


maxima locations was derived. 


4 (C4/C,) sin 8c sin 6 
ae SS ee 

n _ n (VI-1) 
(eeecam F? Forex FT? 


Note that equation (VI-1) is dependent only on the wedge angle 8 and the 


sound speeds of the two fluid media. 


A series of WDGBTM executions were made for varying Co/C, and 8, with 
all other variables held constant. Table (VI-1) is a summary of one case, 
this case is typical of results observed. Results will be discussed in 


Section VI-C. 
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TABLE VI-1l 


Comparison of An from equation (VI-1) and An from WDGBTM results. 


c,/C. = 0.90 


B WDGBTM An An 
maxima from predicted 
locations WDGBTM 

5.0 so nO 7 4.56 
2205 a Oe 3.04 
Seo7, 

nD es 1.40 2733 
55-73 ioe TO 
an. 5 Oot Lose 
Soo 

5 oe alee p29 
3530 0.66 0.88 
3.96 0.49 Oe 
4.45 
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B. MULTIPLE LAMBDA APPROACH 

The second interference pattern calculation also utilizes the 
geometry of Fig. (VI-l). For each c,/C. ratio and beta combination, 
X and As (the wavelength in the bottom media) was calculated according 


to the equations 


a 
= S 
4f sin 6. sing ies 


and 


> 
i 


2 C/E (VI-3) 


where is 150khz. The ratio Ao/KX =O 


is then found and used in the relation 
Z 
vo = [n #(3-x)*]*/? Z [n*+(a-x)7]7/? (vVI-4) 


where N is any positive integer. Note that use of equation (VI-4) 
assumes zero initial phase. 

To find the expected maxima and minima, a plot of 1 vs N was made 
for each 9. Maxima are expected for each integer N, minima for each 
N+1/2. Table (VI-2) summarizes one case of maxima and minima expected 
as a result of equation (VI-4) versus maxima and minima observed on 
WDGBTM output. 

C. SUMMARY OF ANALYSES 
Predictions from both simple analyses are of the same order of 


magnitude as WDGBTM results. The number of maxima and minima are 
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TABLE Via 
OCCURRENCE OF MAXIMA AND MINIMA 


C/C5 = 0.85 


8 (DEG) 1) MAX WDGBTM TY) MAX EXP 1) MIN WDGBTM 1) MIN EXP 
2 4.03 2.200 Dao 4.70 
6.34 
ie Sete Zoe, a5 2 37/5 
OO 4.70 DOG 240 
ere) Zao 2555 2.49 i.7 6 
2aoo b220 S16 
Seow re) 2 Aad? 422) 
B05 4.68 5.40 one 
3.62 5.54 oe D294 
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consistantly predicted as increasing with decreasing wedge angle. 
Occasinal unexplainable phenomena, such as increasing WDGBTM An for 
the 8= 5 on Table (VI-1) and the extra WDOGBTM maxima for 8= 5 
Table (VI-2), were observed, but results are generally supportive of 


an interference effect occurring. 
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VII. COMPARISONS 


A. EXPERIMENTAL DATA 

Netzorg (Reference 4) measured pressure amplitude as a function of 
angle for fixed radii in four media below a wedge shaped upper layer. 
The angle domain version of WDGBTM was used to compare calculated 
pressure patterns with measured patterns. 

All runs utilized finite source input. For each of Netzorg's four 
experimental cases, corresponding computer runs were made with the 
source distance set to 100X, 10X and SX, where S corresponds to the 
actual source distance used in the tank. In all cases S was less than 
10OX. The integrations were carried out over one dumping distance and 
four dumping distances for each source. 

Figure (VII-1) is typical of pressure amplitudes in the bottom. 
The source in the experiment was only about 40 wavelengths from the 
dumping point, even at 10X the source is approximately 110 wavelengths 
away. Since the method of images calculations utilized to determine 
the pressure amplitude along the wedge bottom assume a plane wave, it 
is questionable that this program can produce a realistic comparison 
with the measured data. 

Figures (VII-2) and (VII-3) illustrate the wide variation in results 
obtained depending on the source distance and integration interval. 
Calculations demonstrate expected trends, l.e., pressure amplitude 


increases with decreasing radius distance (9), with closing source 
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distance, and as a function of an increasing integration interval over 
the wedge bottom. The beam depression angle becomes more shallow as 
source distance increases. Results at very low angles are not 
considered valid because of assumptions inherent in the WDGBTM model. 

The circled data points of Figures (VII-2) and (VII-3) are a 
transposition of data from Reference 4. Since absolute amplitude was 
not specified in the reference, rescaling to the case most near to the 
actual conditions was performed, i.e., the data ee normalized to the 
case x, = 1.0 and srcd = 3.5X. The shape of the experimental curve and 
the apparent beam depression angle are the only valid comparisons that 
should be made. 

Comparisons of WDGBTM beam depression angles with measured values 
gave varying results. Measured depression angles decreased with 
increasing 0, WDGBIM depression angles increased. WDGBTM results were 
within 6 percent to 62 percent of measured values. No correlation was 
observed between prediction accuracy and source distance or prediction 
accuracy and the interface integration interval. All measured beam 
depression angles in Reference 8 were for the two shortest radii of 
each case, the three longer radii could only be qualitatively compared 


with graphical data Possible attenuation effects were not considered. 


Numerical results are contained in Appendix D. 


B. PARABOLIC EQUATION APPROACH 
Kuperman (Reference 9) describes a Parabolic Equation model result 
for a wedge problem, and compares these results to one set of data 


extracted from Reference 10. 
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The Parabolic Equation results show very well defined beams in the 
bottom, with no apparent interference effects as observed in comparable 
WDGBTM runs (Figures V-1 and V-4). kKuperman used an attenuation factor 
of 0.5¢dB per wave length, a value equivalent to an input of 0.984 
Nepers/X for WDGBTM. It should be noted in Figure V-4 that there is 
complete separation of the beams for WDGBTM when &= 1.0. 

Kuperman states agreement within 20 percent of a measured data set 
of Reference 10. This measured data is identical to Case 1 of Appendix 
D. WDGBTM accuracy for this case varied from zero percent error to 22 


Bemcent error with error calculated according to the formula 


(observed data) - (measured data) 


100 x Measured data 
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Vir. CenecEusrtons 


1. Program WDGBTM offers a fast, efficient means of analysing the 
behavior of sound in a fast bottom underlying a wedge-shaped medium. 
Effects of attenuation are easily demonstrated, and modes desired are 
readily controlled by specification of source location and integration 
interval selected. Results of the model are not inconsistant with 
current measured data, and illuminate observations generated by the 


Parabolic Equation model. 


2. A simple equation for calculating beam depression angle as a 
function of wedge characteristics was obtained which offers sufficient 
accuracy for use as a rapid analytical tool. The equation was derived 
for acoustic sources placed at one-half the wedge angle. Accuracy degrades 
when the wedge medium density to bottom medium density ratio decreases 
On U.S. 

3. It is recommended that further experiments be performed that 
would further test the model. New experiments should be sufficiently 
large scale to ensure plane waves at the wedge bottom interface and to satis- 
fy the k,|r-r° |>>2n eri teria. Bottom media that reasonably model 


real ocean fluid-like bottom are also recommended. 
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